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We estimate the growth of matter perturbations in a class of recently proposed dark-energy 
models based on the (loop-corrected) gravi-dilaton string effective action, and characterized by a 
global attractor epoch in which dark-matter and dark-energy density scale with the same effective 
equation of state. Unlike most dark-energy models, we find that the accelerated phase might start 
even at redshifts as high as 2: ~ 5 (thus relaxing the coincidence problem), while still producing at 
present an acceptable level of matter fluctuations. We also show that such an early acceleration is 
not in conflict with the recently discovered supernova SN1997ff aX z ~ 1.7. The comparison of the 
predicted value of erg with the observational data provides interesting constraints on the fundamental 
parameters of the given model of dilaton-dark matter interactions. 



I. INTRODUCTION 



Independent cosmological observations have recently pointed out the existence of a signiflcant fraction of 
critical density in the form of unclustered matter, possibly characterized by a negative pressure 0, ^, ^, ^. 
Such a component of the cosmic fluid, conventionally denoted dark energy, is believed to drive the accelerated 
evolution of our present Universe, as suggested by the study of the Supernovae la (SNIa) Hubble diagram. 
While the simplest explanation of the dark-energy component is probably a cosmological constant, it might 
be neither the most motivated, nor the best fit to the data. 

The existence of a (presently dominating) dark-energy component raises at least three important questions: 
a) what is the physical origin of such a new component? h) why its energy density today is just comparable 
to the matter energy density? c) did the acceleration start only very recently in the cosmological history? 

The first question refers to the fundamental physical mechanism able to generate a cosmic and universal 
dark-energy distribution. So far, most models of dark energy have adopted a purely phenomenological 
approach, even if a few proposal concerning the_possible role of the dark-energy field in the context of 
fundamental physics have already appeared P, |^, 0, |, ^ [ict] . 

The other questions concern two distinct (although possibly related) aspects of the so-called "coincidence 
problem" jl^]. The first aspect refers to the density of the dark-energy and of the dark-matter component. 
Since in most models the two components have different equations of state, they scale with time in a different 
way, and their densities should be widely different for essentially all the cosmological history: by contrast, 
current data are strongly pointing at a comparable proportion of the two components just at the present 
epoch. 

The second aspect of the coincidence problem may arise if the epoch of accelerated evolution started 
only very recently, say &t z 'p^ 1. There are reasons to suspect that this may be the case because, in most 
models, the growth of structures is forbidden after the beginning of the accelerated expansion, so that the 
acceleration cannot be extended too far into the past. It has also been recently claimed that the supernova 
SN1997ff ||l^, |l^ at 2: ^ 1.7 provides an additional indication that the acceleration is a relatively recent 
phenomenon. 

A promising approach for a possible simultaneous answer to the above questions has been recently provided 
by a dilatonic interpretation of the dark energy based on the infinite bare-coupling limit of the superstring 
effective action Q, whose cosmological solutions are characterized by a late-time global attractor where dark- 
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matter and (dilatonic) dark-energy densities have an identical scaling with time. A very similar cosmological 
scenario has been studied also in |0, |l^. The respective amount of dark-matter and dark-energy density 
is eventually determined by the fundamental constants of the model, and it is expected to be of the same 
order, so that their ratio keeps frozen around unity for the whole duration of such an asymptotic regime. 
In such a context, the today approximate equality of dark-matter and dark-energy density would be no 
longer a coincidence of the present epoch, but a consequence of the fact that our Universe has already 
entered the asymptotic regime (actually, for a region of parameter space allowed by observations, it is also 
possible that the ratio of dark-energy to dark-matter density is close to unity not only in the asymptotic 
regime, but already after the equivalence epoch). Finally, the cosmic field responsible for the observed large- 
scale acceleration would be no longer introduced ad hoc, being identified with a fundamental ingredient of 
superstring/M-theory models of high-energy physics. 

In this paper we will focus our discussion on the possible constraints imposed by structure formation on 
the above model of dilatonic dark energy, by studying the growth of matter perturbations in the two relevant 
post-equivalence epochs: a first decelerated epoch, dubbed dragging phase, and a second accelerated epoch, 
dubbed freezing phase. We will then reconstruct the behavior of the cosmological gravitational potential, 
evaluate the Sachs- Wolfe and integrated Sachs- Wolfe contributions to the spectrum of CMB anisotropics, 
and estimate the present level of matter fiuctuations (in particular, the so-called variance as, smoothed out 
over spheres of radius 8 Mpc h~^). Finally, we will compare it with observations. 

The results that we find are interesting, and might help answering the third question posed at the beginning 
of this section: among the range of parameters compatible with a phenomenologically acceptable erg we find 
indeed values allowing a long epoch of acceleration, starting as far in the past as at z ~ 5. By contrast, the 
models of dark energy uncoupled to dark matter, with frozen (or slowly varying) equation of state, cannot 
accelerate before z ~ 1 . It seems appropriate to anticipate here that the production of an acceptable level of 
fiuctuations even in the case of an early start of the accelerated epoch is due to two concurrent factors: the 
first is that, during the freezing phase, perturbations do not stop growing like in other models of accelerating 
dark energy; the second is that the horizon at equivalence in our model shifts at larger scales with respect 
to a standard ACDM model. It is also important to stress that such an early acceleration is by no means in 
contrast with the recently observed supernova SN1997ff at z « 1.7. 

For the model of dilatonic dark energy considered in this paper the coincidence problem can thus be 
alleviated by the fact that the ratio of dark-matter to dark-energy density is of order one not only at 
present, not only in the course of the future evolution but (at least in principle) also in the past, long before 
the present epoch. It should be clear, however, that this possibility could be strongly constrained by future 
observations of supernovae at high redshift. 

The paper is organized as follows. In Sect. II we present the details of our late-time, dilaton-driven 
cosmological scenario, and define the (theoretical and phenomenological) parameters relevant to our com- 
putation. In Sect. ^ we discuss the growth of matter perturbations and compute the variance as- In Sect. 
uM we impose the observational constraints on our set of parameters, and determine the maximal possible 
extension towards the past of the phase of accelerated evolution. We also find, as a byproduct of our analysis, 
interesting experimental constraints on the fundamental parameters of the string effective action used for 
a dilatonic interpretation of the dark-energy field. In Sect. |^ we compare our model with the constraints 



provided by the farthest observed supernova at z ~ 1.7. Our conclusions are finally summarized in Sect. VI 



II. THE MODEL 



The model we consider is based on the gravi-dilaton string effective action ||l6| which, to lowest order in 
the higher-derivative expansion, but including dilaton-dependent loop corrections and a non-perturbative 
potential, can be written in the String-frame as follows 1^ 



-4,(4,) 



R + Z 



+ Sm{g, 4>, matter) 



(1) 



Here Xs is the fundamental string-length parameter, and the tilde refers to the String-frame metric. The 
functions ip{(j)),Z{(l)) are appropriate "form factors" encoding the dilatonic loop corrections and reducing, in 
the weak coupling limit, to the well known tree-level expressions Jl^ ] Z = exp(— ?/j) = exp(— (^). 

In this paper we are interested in a "late-time", post-big bang cosmological scenario, in which the dilaton is 
free to run to infinity rolling down an exponentially suppressed potential, and the loop form factors approach 
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a finite limit as 4> ^ +00. Assuming the validity of an asymptotic Taylor expansion and following the 
spirit of "induced-gravity" models in which the gravitational and gauge couplings saturate at small values 
because of the large number (A*' ~ 10^) of fundamental GUT gauge bosons entering the loop corrections, we 
can write, for (/> — > -|-cx3, 

g-'/'W ^ cj + bie-'^ + 0{e-^'^), Z{4>) = -cl + h2e-'* + 0{e-^'>') , 

V{4>) ^ Voe-'" + 0{e-'"l'). (2) 

The dimensionless coefficients cf,c| are typically of order lO'^, because of their quantum-loop origin. We 
may note, in particular, that cf asymptotically controls the fundamental ratio between the (dimensionally 
reduced) string and Planck scales 0, cf = A^/Ap, which is indeed expected to be in the range 10 \p/\a — 
0.3-0.03. 

To complete the model we have to specify the matter action Sm of eq. (|l]), containing the coupling 
(possibly renormalized by loop corrections) of the matter fields to the dilaton. The variation of Sm with 
respect to (j) defines the (String-frame) dilatonic charge density ct, whose appearance is a peculiar string 
theory effect p9| |, and represents the crucial difference from conventional (Brans-Dicke) scalar-tensor models 
of gravity. 

For the cosmological scenario of this paper we shall assume, as in |^, that Sm contains radiation, baryons 
and cold dark matter, and that the dilatonic charge of the dark matter component switches on at sufficiently 
large couplings, being proportional (through a time-dependent factor q{4>)) to its energy density pc- Also, 
q{(j>) is assumed to approach a constant (positive) value go as ^ -l-oo, 

g(0) = go-fO(e'^«^). (3) 

The dilatonic charge of radiation and of ordinary baryonic matter are instead exponentially suppressed in the 
strong coupling regime Q, and this guarantees the absence of unacceptably large corrections to macroscopic 
gravity, since in the model we are considering the dilaton is asymptotically massless and leads to long-range 
scalar interactions (see however |^ for possible testable violations of the equivalence principle, and other 
non-standard effects, possibly observable in such a context). 

By considering an isotropic, spatially fiat metric background, and a perfect fiuid model of matter sources, 
it is now convenient to write the cosmological equations for the action (Q) directly in the Einstein frame, 
defined by the conformal transformation g^^ — cf exp(i/'). In the cosmic-time gauge (and in units 
2Ap = l&nG — 1) the equations are |^ 

^ = p+p^ , 4H + 6H^ = -p-p^, (4) 
fe'(0) {4> + 2.Hi) + k{<f>) k'{^) <j>'' + V'{<t)) + i [i''{^){pt + Pc) + q^Pc] = 0, (5) 

where a prime denotes differentiation with respect to (p, and 

fc'(0) =3^''-2e2'^Z, V = cie^^V, q{4>) ^ a, / p,, (6) 
P4>=\k\4.)4>^ + V{(p), p^^'^k'{<t>)4? -V{4>), (7) 

p=Pr+Pb+Pc, p = Pr/3. (8) 

We have explicitly separated the radiation, baryon and cold dark matter components {pr, pt, Pc), and intro- 
duced the (Einstein-frame) dilatonic charge per unit of gravitational mass, q{4>), which is non- vanishing (at 
large enough 0) only for the dark-matter component. The combination of the above equations leads to the 
separate energy conservation equations: 

Pr + 4:Hpr = 0, (9) 
Pb + 3Hpb - ^(ptp'pb = 0, (10) 

Pc + 3Hpc - + q)pc = 0, (11) 

p^ + 3H{p^+p^) + ]-4> [^'{pb + Pc) + qPc] = 0. (12) 
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With the above assumptions on tp, Z and q, and for appropriate values of the parameters of the loop 
functions (in particular, for a sufficiently small value of Vo), it has been shown in |^ that the phase of 
standard, matter-dominated evolution is modified by the non-minimal, direct coupling of dark matter and 
dilatonic dark energy. After the equivalence epoch, in particular, the Universe may enter a phase of "drag- 
ging", followed by an accelerated phase of asymptotic "freezing". For the reader's convenience we recall here 
the main properties of such two phases, referring to j^, Q for a more detailed discussion (see also |Q for a 
general study of the dynamical system). 

(1) "Dragging" phase . The potential V{(j}) is negligible, the evolution is decelerated, P0 is still subdominant 
(as well as pb), but evolves in time like the dominant component pc, so that the dilaton dark energy is 
"dragged" along by the dark matter density. 

In this phase k{(p) can be approximated by a constant, k{(f>) = fci, and it is thus convenient to define the 
rescaled field (p — ki(j) which has canonical kinetic term in the action, and which satisfies with Pc the system 
of coupled equations 

1 f tp' + q 



' + |^^-^jp, = 0, (13) 

1 f i>' + q 



We can then define the canonical effective coupling of the dilaton to dark matter by the function /3(0), 
defined by: 

which, in the dragging phase, is also approximated by a constant, /3((?!)) = /3i <C 1 (the conventional factor 
l/\/3 has been introduced here to adapt the notations of this paper to previous studies of the dark-matter- 
scalar system |2l||). 

Using eqs. (J), ( [l^ ) we find, in the dragging phase |^, 

^=-4Hl3i/V3, (16) 

so that, from eq. ([T^, 

p,^i/2^p,^a-(^+^''?/3). (17) 

Because of the dragging the time evolution of the dark matter density deviates from the typical behaviour 
of dust sources, in such a way that pc decays slightly faster than energy density of baryons, pb ~ a~^. It is 
however unlikely that this effect may lead the Universe to a baryon-dominated phase, because this trend is 
soon inverted in the subsequent, freezing phase. 

(2) "Freezing" phase . The asymptotic dilaton potential V{(j>) — Vo exp(— (^) comes into play, the evolution 
(for large enough values of g) is accelerated, the dark matter, potential and (dilatonic) kinetic energy densities 
evolve in the same way, so that the ratio p^/pc is frozen to an arbitrary constant value. The critical fraction 
of dark matter, potential and kinetic energy densities are also separately constant throughout this phase 

In such a phase, asymptotically approached when (p +oo, one has q{(l>) = go, tp' = 0, and the parameters 
k((p), f3{(p) can be again approximated by constant values k2, 13-2 (in general different from the previous ones), 
related by fe2/32 = V3qo/2. The coupled equations for the canonically rescaled field (p = k2(p are modified by 
the presence of the potential 

^ + SH(P+^ + ^Pc = 0, (18) 
d(P V3 

Pc + 3Hpc- ■^Pc'$=0, (19) 
v3 

and, together with eq. (^, are solved by the following configuration |^: 

V = Voe-^/"^ ^ fl-'^/^^+W) . (20) 
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In units of critical energy density we find, in particular, 
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Note that, from eq. (po[). 



2 + go' 12//2 (2 + go)2' 

n0 = f7v + f^fc, He = 1-^0. (21) 

so that the freezing phase is accelerated for go > 1- In that case, according to eq. ([20[), the dark matter 
density tends to be strongly enhanced (as time goes on) with respect to the baryon density, which on the 
contrary is uncoupled to the dilaton and thus evolves in the standard way, pb ~ a~^. As already pointed 
out in ig, [i^], it is tempting to speculate, in such a context, that the present smallness of the ratio pb/pc 
could then emerge as an artifact of a long enough freezing phase, started before the present epoch. 

In this paper, in order to discuss the possible bounds imposed by present observations on the above 
scenario, we will consider a simplified model of late-time (i.e., after-equivalence) cosmology, consisting of 
two phases. More precisely, we will drastically approximate the background evolution by assuming that the 
Universe performs a sudden transition from the radiation-dominated to the dragging phase at the equivalence 
epoch a = aeq, and from the dragging to the freezing phase at the transition epoch a = af. We will discuss 
in this context the phenomenological constraints on the parameters of the string effective action, and in 
particular their possible consistency with an early beginning of the freezing epoch, zf = (ao/af) — 1^1, 
which (as already mentioned) may be relevant for a truly satisfactory solution of the coincidence problem. 

To make contact with previous results we will use here the following explicit model of dilaton potential, 
charge and loop corrections (already adopted in j^): 

g-V>(« ^ + ,2^ z{^) = e-^ - cl q{<t>) = go T, (23) 

V{(l)) = Cl mve^''' [exp (— e~'^/Qi) — exp (— e~'*/Qf2)] , < Q2 < ai. (24) 

Here Ci,C2, c? are numbers of order 10^, the asymptotic charge satisfies go > 1 to guarantee a final accelerated 
regime, and the dilaton potential reduces asymptotically to the exponential form of eq. (^, with Vo = 
my{oii — a2)/aia2 (in units 2Ap — 1). It is worth noting that the mass scale mv , controlling asymptotically 
the amplitude of the potential (and thus the beginning of the freezing phase), is possibly expected to be of 
non-perturbative origin, and thus related to the fundamental string scale in a typical instantonic way. 



mv = exp 



Ms, (25) 



where Q^quT — ^/^^ '® asymptotic value of the GUT gauge coupling, and (5* is some model-dependent 
loop coefficient. As noted in ||9|, a value of /3* slightly smaller than the usual coefficient of the QCD 
beta-function is already enough to move mv from the QCD scale down to the scale relevant for a realistic 
scenario of dark-energy domination; a typical reference value is, for instance, /3* ~ 0.36, which corresponds 
to mv ~ Ha ~ 10"^^ Mp, and thus to a freezing phase starting around the present epoch. 

With the above explicit forms of the loop corrections we can now compute the constant parameters for 
the dragging (fci,/3i) and the freezing (fc2,/32) eras. By setting /i2 = \/3ci/(\/2c2) we find 

, \/3 g0M2 u R 90^2 (0R\ 

ki = , Pi = -TTT' ^2 = , P2 = (26) 

P2 2c2 fl2 2 

The constant fi2 / Vs represents the slope of the dilaton potential V{(l>) during freezing. With these definitions, 

go = 2/32//i2,c' =/32//?l. 

For our model of background we can finally express the phenomenological variables, required for the 
subsequent computations, in terms of the above set of parameters. From eqs. ([l7[), ( po| ) we get the barotropic 
parameter w = {ptot/ptot) + 1 relative to the effective equation of state of cold dark matter. 



(27) 
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in the dragging and freezing eras, respectively. Another useful parameter is ilc- in the dragging phase, from 
eq. (0), 



In the freezing phase, from eq. 



12m 



1 



+ 2/32 ~ 9 



(28) 



(29) 



2(/32+M2)2 

(this is therefore to be identified to the value observed today). 

Note that the ratio of dark-energy to dark-matter density may be close to unity throughout the cosmic 
evolution after equivalence, if the constants Pi, [32,^.2 are of the order of unity. In this sense, in such a model, 
it is even possible that a serious coincidence problem never really arises. By contrast, dark-energy models 
without a stationary phase in which pm ~ (see e.g. |^, |24[) have to explain the today value (or order 
one) of ratios Q,^ /i^m which range from extremely small values in the past to unboundedly large values in 
the future Q. 

We end this section with the computation of three useful red-shift parameters ZE,Zf,Zb corresponding, 
respectively, to the radiation-matter equivalence, to the beginning of the freezing epoch, and to the baryon 
epoch (associated to a possible baryon-dominated phase). Prom the behaviour of and pt in the freezing 
phase, 



.(ao) (^) 



3W2 



pb = Pi,(ao) 



P0 = P,p[ 

we can determine the baryon red-shift epoch a;,, such that p,/, = pt, as follows |22 



at \ilt 



3/32 



(30) 



(31) 



The present ratio (f2^/r2i,)o is a known observational input. Our model excludes the possibility of a baryon- 
dominated phase and thus requires, for consistency, zi, > Zf (see p2| for a dark-energy model with a baryonic 
epoch). 

The equivalence scale can be obtained by rescaling pc and pr from down to ao, i.e. 



Pc(ao) = , 









(-) 







from which 



1 I 

1 + Ze = = 

tte 



3{w2—wi) 



(32) 



(33) 



where, again, (r2r/f^c)o is an observational input. Similarly, we can determine the freezing epoch by rescaling 
the dilaton potential energy, 

^ ao 



pv{af) = Voe = pv(ao) 



(34) 



Here pv{ao) = GH^Qv, where f2v is determined by eq. (Ell), and cj>f is the transition scale between small 
and large values of the dilaton charge, namely (j)f = (2/go)Tnc, from eq. (Esl). We thus obtain 



1 I "0 

1 + Zf = — = 

a/ 



Vo 



M2 +/^2 
3m2 



9 + 4(32{p2 + (32) 



(35) 



(36) 



4(p2-H/32)2 ' 

where Ho is the value of the Hubble parameter provided by present observations. 

In conclusion, we have a four-parameter model of background, spanned by two possible (related) sets 
of independent variables: the phenomenological variables {l3i, l32, P2, Zf} or, equivalently, the fundamental 
variables {ci/c2,c^ ,qo,Vo}, referred to the (j> — > -|-cxd limit of the string effective action of eq. (|l|). The 
mapping between the two sets is defined by eqs. (Ea) and (pd). Note that two of the four parameters can 
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be in principle determined by fitting the observational data relative to the present fraction of cold dark 
matter, iln, and the dark-energy equation of state, W2: we can determine, for instance, /i2 and P2 through 
eqs. (Eq), (pU)- In the following sections we will discuss the allowed regions left by various phenomenological 
constraints in such a parameter space. 



III. CONSTRAINTS FROM STRUCTURE FORMATION 



In this Section we will compute the r.m.s. dark-matter density contrast erg for the model presented in 
Sect. II, combining it with the CMB angular power spectrum (at low multipoles) in order to eliminate the 
dependence on the normalization factor. Another analytical computation of as including dark energy has 
recently been presented in ||2^ ], but only for models in which dark energy and dark matter are uncoupled. 

Before starting the computation, let us recall some preliminary condition to be imposed (for consistency) 
on our model of background. 



A. Consistency conditions on the background 



As illustrated in Sect. II, the evolution of our cosmological background is characterized by two stationary 
{pc ~ P4,) stages: the dragging era (labelled by the subscript 1), and the final accelerated freezing era 
(labelled by 2). This scenario can be consistently implemented provided the parameters are chosen in such 
a way as to satisfy the following conditions. 

First, the dragging era exists, and is a saddle point of our dynamical system, only if 

< \/3/2 ~ 0.87. (37) 

Second, in order for the baryons not to dominate (in the past) over the dark-matter component, it is necessary 
to impose (as anticipated): 

zt > Zf. (38) 
Third, the final freezing phase is accelerated only if: 

At2< 2/32. (39) 

The current SNIa observations, however, require more than simply an acceleration, as shown by a recent 
analysis |^ of SNIa data in models which include a freezing epoch. The result is that only models with 
an effective equation of state W2 ~ 0.4 (best fit) or W2 < 0.5 (at one sigma) are consistent with the SNIa 
Hubble diagram. 

In the following discussion we will reduce, for simplicity, our set of free parameters by using the observed 
CDM fraction of critical density. Sic — 0-3, as a fixed observational input. This can be used to eliminate fj,2 
through eq. (p9|), in such a way that we are left with three parameters only, /3i,/32 and Zf. In that case, the 
above limits on W2 = ^2/(^2 + ^2) define two reference values for /32, 

P2 = 4.02, (best fit ) (40) 
P2 = 2.35, (one sigma ), (41) 

which will be used in our subsequent analysis. For future reference, note that the lower limit at 95% c.l. 
is f32 > 1.55. We note, finally, that pushing back in time the transition between the dragging and freezing 
eras, that is increasing Zf, implies a decrease of Ze- The condition 

Zf < 100 (42) 



prevents the unwanted crossing of these two quantities. 

All the above constraints will be imposed on our subsequent computations. 
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B. Angular power spectrum at low multipoles 



In this subsection we will extend the treatment of the Sachs- Wolfe effect presented in |27[ (to which we 
refer for the notation) in such a way as to include the case of coupled dark energy and dark matter. We shall 
assume a conformally flat metric, ds^ = a?{dT^ — dx^dxi), and we shall exploit the fact that, in the absence 
of anisotropic stress, the two scalar potentials ^' and <& defined in the longitudinal gauge turn out to be equal 
and to coincide, in the Newtonian limit, with the usual gravitational potential. The general expressions for 
the Sachs- Wolfe (SW) and integrated Sachs- Wolfe (ISW) parts of the angular power spectrum, for adiabatic 
scalar perturbations, can then be written as [p^ : 



= ^ I dkk^\-f{k,ra)\^jj[k{ro-Ta)], (43) 
= 2Gi(£) / dk\<f'{k,n)\\ (44) 



where G2(£) = r[(£ + l)/2]/r[{£ + 2)/2], ji are the spherical Bessel functions, the prime denotes differentia- 
tion with respect to conformal time, and Td is the conformal time at decoupling, while n — tq — {i+ l/2)/k 
(notice that the above expression for the ISW coefficient has been already integrated over r). 

Let us denote with 5k the CDM density contrast for the wavemode k, and with 4> and x = ^4') respectively, 
the values of the background scalar field and of its fiuctuation. The label will denote the present time to, 
and it is to be understood that all the density parameters ^lc,^b, etc. are always referred to the present 
time, unless otherwise stated. Finally, all perturbation variables will be expressed in terms of their Fourier 
components. The potential ^ is then determined by the relativistic Poisson equation (in units IGtvG — 1) 

*W - ~i,H!n.Mao)sik,a) _ ^x'? - ^X? - ^x'^, (45) 

where we have introduced the (possibly fc-dependent) function s{k,a), which accounts for the time-evolution 
of the dark matter and of its density contrast, according to the equations 

Pc{a) = pc{ao)f{a), (46) 
5fe(a) = 5k{ao)D{k,a), (47) 

3{k,a) = D{k,a)f{a) ( — V. (48) 
\ao/ 

The quantities appearing in the r.h.s. of eq. (jlSj ) are evaluated in the synchronous gauge, while the 
potential ^' is calculated in the longitudinal gauge since this helps extending the validity of the Poisson 
equation to all scales. Note that we have dropped from the Poisson equation the contribution of the dark- 
matter velocity fiuctuations, since it can be shown that they are negligible ||l5|| . The baryon contribution 
has been neglected as well. Defining the matter power spectrum P{k) = |(5fc(ao)|^, it follows that 



Cf^ = ^i/o*nc(l + S) H^P{k)\s{k,adtjj[k{To-Td)l (49) 



Cl''' = iH^^lGlmi + Il) J ^P{k)\s'[k,a{r,)]\ 



(50) 



The functions E and 11 represent corrections due to the scalar field contribution. In the next few paragraphs 
we will neglect these effects, concentrating the attention only on the dark-matter part of the power spectrum, 
while the scalar contribution will be reconsidered at the end of our calculation. The convenience of this 
procedure will become apparent later on. 

Focusing our attention on the background evolution after the epoch of matter-radiation equivalence, we 
can now specify the evolution of the CMB energy density by setting 

/(a) = {a/ao)-''"\ a>af (51) 
= {af/ao)-'^"'{a/af)-'-\ a<af. (52) 



where wi^2 were defined in eqs. (|27| ). For what concerns the growth of matter perturbations, since we are 
interested in the low-multipole branch of the spectrum, it is reasonable to consider only scales that reenter 
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the horizon after equivalence (and before freezing, as in the accelerated stage no reenter is possible). The 
growth of Sk has been derived in |15[ as a function of the parameters /3 and fi of the dragging and freezing 
epochs. It turns out that the evolution is the same for all modes: 



D(k,a) = (a/ao)™^, a>af 

= (a/a/)'"i(a//ao)'"', a < a/, 

where 

1^4 2 A - 10/32 - 

= m2 = ^^^-^^, (53) 

and 

= -108 + 44/32M2 + S2|3|^i2 + 25At2 + Pi (32^2 - 44) (54) 

(in the dragging phase the same exponent mi describes the growth of perturbations both inside and outside 
the horizon). Considering scales with k < kd, where the subscript d stands for decoupling, the relevant 
function for the computation of the ordinary SW effect can thus be parametrized by a fc-independent function 
as follows: 

s(a/) = (ad/a/)'"i~^"'i(a//ao)'"^"^"'^(ad/ao)^ = (af/ao)"^. (55) 
where 0:2 = rn2 — 3w2 + 2 (note that Qi = mi — 3wi +2 = 0). 

Some comments are now in order, concerning the evolution of perturbations described by the above 
equations. First we notice that in the dragging phase mi > 1, i.e. that the perturbations growth is faster 
than in a standard CDM model. This is due to an extra pull on dark matter arising from the dark-energy 
coupling, that act as an additional (scalar) gravity force Secondly, the growth in the accelerated 

freezing regime does not vanish asymptotically, like in other dark-energy models. This is again an effect 
of the dark energy-dark matter coupling and of the fact that the dark-matter density is not driven to zero 
by the acceleration. Finally, since the baryons are decoupled from the dilaton, they evolve differently, and 
a bias between the baryons and the dark-matter distribution is expected. The constraints from this effect 
have been discussed in ref. jisf . 

It is important to stress that the matter power spectrum, calculated today for scales k < ke, where fee is 
the scale that reenter the horizon at equivalence, 

ke ^ a^He ^ aeHo [(a//ao)~^™^(ae/a/)~^™^] ^''^ , (56) 

does not change respect to the primordial shape 

P{k) = Ak", (57) 

where A is the usual normalization factor (see e.g. ref. |^^ ). This is due to a peculiarity of the dragging 
phase, i.e. to the fact that the perturbation growth is identical inside and outside the horizon, just like 
in the case of standard CDM models (although, as already mentioned, the growth rate mi is larger than 
unity). The dragging phase, in addition, leads to a time-independent gravitational potential (again, like in 
standard CDM models), and thus it does not contribute to the ISW effect, which is entirely produced in the 
subsequent freezing phase. 

Prom the integration of eq. (^9|), with E = we then easily obtain 

C!"^ = ^AHS+'nls\af)G{£,n), (58) 
Id 

where 

GU n) - r(3-n)r((2£ + n-l)/2) 

and where the factor 7 appears by eliminating ro (from the result of the SW integral) in terms of the present 
Hubble scale -ffo, according to the definition 

f da 2 , ^ 
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This integral can be finally estimated by considering the separate contributions from the two phases of our 
model, and we obtain: 



2X2 



^ 3 Aw jaf/ao) 



2X1 



3 3 

Ai = -Wi — 1, A2 = -W2 — 1, Aw = W2 — Wl 



(61) 



(the usual result for the standard cosmological model is instead 7 = 1). 

It is appropriate to reconsider at this point the scalar-field corrections to the SW effect, represented by E. 
The scalar field fiuctuations which are outside the horizon during the dragging phase grow proportionally to 
the CDM density contrast (see ||2^, |2|]). Since pc ~ P0, it is found that E is a constant, and depends only 
on /3i as follows: 



6VfA^(3 + /3f)(15 + 4/30 

405 + 252/32 - 336/3* + 64/3f ^ ' 

(the contribution coming from the scalar-field potential has been neglected, since in the dragging phase the 
dilaton kinetic energy is dominant with respect to V). 

For what concerns the ISW effect, the scalar field contribution on sub-horizon scales is negligible | p2[ , 
and we can consistently set II = 0. If we define the variable y — fcro, and we use the result that, in the 
accelerated epoch, 

s(a) = (a/ao)"^ (63) 

then we are lead to 

9 



The lower limit of integration, 



i — r To ^7Ai 



(64) 



(65) 



has been obtained by imposing T( > r/, since only in the second (freezing) stage there is a significant ISW 
contribution. The time-dependence of the freezing scale-factor, on the other hand, can be parametrized by 

a = ao [X2 (Hor + B)]^^^' , 

where 



(66) 



\aoJ VA2 Ai 



(67) 



By using all the above results the ISW integral can be finally performed and, in the case n = 1, the result 



r 9 [(i? + 27)A2]^^ + [(i? + 27r)A2]^^2 K [4^a^(r - 1) - {B + 2'yr)X2] 
C,,jsw = -An^G2moa2 . (68) 

For a generic, primordial spectral index with n 7^ 1 a much more complicated (but still analytic) expression 
may be obtained. 

In conclusion, the dimensionless, angular power-spectrum at low multipoles {£ < 10) can be approximated 

by 

Cf = Ce.sw + Cejsw = AF, (69) 

where, for n = 1, 



F{£;af, (3,, 132, fi2)=Ht,n', 



+ ^Gl{e)a2 



47v£{l+l) (S) ^^^^^ 

[{B + 27)A2]'^ + [{B + 27r)A2]''^~' ■ [47Q2(r - 1) - (B + 27r)A2] 
72(2^ 1)2 (2a2-A2) 



(70) 



The variable Ci , representing the experimentally observed angular power spectrum, measured in units of 
{iiKY , can thus be finally written in the form 

TSi{£+l)AF^ ^^^^ 



2n 



where To = 2.726 x lOVK. 
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C. Calculation of ag 



The (dimensionless) variance of the CMB density fluctuations, in spheres of radius Rs = 8h ^ Mpc, where 
h = //o/(100kmsec-^Mpc-^), is defined by ||| 

"i^^ j Pik)Wiik)k'dk = A.R-(^+^)i,^ (72) 

where Ws (k) is the spherical top-hat window function of radius -Rs , and 

Ii= [ x'+"T\x)W^(x)dx. (73) 



Note that we are using the full power spectrum corrected by the transfer function T{k), i.e. P{k) — 
Ak"'T^{k), since in the definition of as it is necessary to include also scales with k > fee. 

It is important to stress, at this point, that the above transfer function is identical to the one of the 
usual ACDM model. During the dragging phase, in fact, the perturbation growth does not depend on the 
wavenumber, while during the freezing phase only sub- horizon perturbations have to be taken into account, so 
that no distortion of the power spectrum occurs after the equivalence epoch. Therefore, the transfer function 
only expresses the usual correction to the primordial spectrum due to the different growth of perturbations 
in the radiation epoch (those entering the horizon before equivalence are depressed with respect to those 
entering later). Since, in our model, the cosmological evolution before equivalence is standard, we can 
safely adopt the transfer function of a ACDM model for which the wavenumber ke crossing the horizon at 
equivalence is the same as in our model. For a ACDM model with present density ^.^.(^/^c^jyi-) one has, in 
particular, 

ke{Acdm) = 0'OHoQc{Acdm.)y (74) 

By equating to fee(Acdm) the value of fee determined in our model (see eq. we thus obtain the effective 
density parameter 



to be used for the determination of the equivalent transfer function (we shall of course restrict our analysis 
to the case iic(Acdm) < !)• For the final numerical integration of eq. ( {73| ) we shall use the ACDM transfer 
function proposed in I pof . 

Combining eq. ( [ti] ) and eq. ([t^), in order to eliminate the normalization factor A, we finally obtain 



2 _ Cfll , , 

" 7rr2i?|£(^ + l)F' ^ ' 

where Ii depends on OajAcdm.)- The observed C( has been obtained by fitting the COBE data as in [ po| . 
The comparison of eq. ( |7q ) with the experimental value of ag, 

(78 = (0.56 ±0.1)^7J°■''^ (77) 

taken from data on clusters abundances will eventually give the sought-for constraint on the parameters 
of the dilaton model introduced in Sect. 11. 



IV. RESULTS 



In order to implement the constraints imposed by as we shall first eliminate /i2 from our set of parameters 
(as already anticipated), by using eg. ( |29| ) and fixing the CDM density at the value fie = 0.3, suggested by 
present observation (see e. g . P2]). Also, we shall assume for the baryon density the value predicted by 
standard nucleosynthesis |33| , Qb = 0.02 h~^, and we set h — 0.65. 

The value of (32, at given /i2, should be unambiguously determined by go, and then by the observed 
value of the cosmic acceleration through eqs. (|2^) and (|2^). However, in view of the present experimental 
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uncertainties, we have accepted here an open range of possibilities and we have illustrated the constraints at 
two (rather different) values of fi2- ~ 2.35 (the minimum value allowed at one sigma by SNIa and 
132 ~ 4.02 (the best fit to the SNIa data |^^)). These two reference values will be used in all the following 
discussion. 

We shall first illustrate the erg constraint in Fig. 1 by plotting the curves corresponding to the experimental 
values ( [ttI ) (with a 3 sigma error band) in the plane {zf, f3i} with (32 fixed . The upper curves (and the 
darker regions) corresponds to lower values of as- In this figure (and in the next one) the white region has 
been excluded because the transfer function parameter i^c(Acdm) becomes larger than unity. The dashed 
horizontal lines represent the upper bound on 2:/ imposed by the baryon constraint Jj^). The allowed region 
is below the dashed line, and within the upper and lower white curves. In the left panel of Fig. 1 we have 
used /32 = 2.35, and we obtain that the maximum past-extension of the accelerated regime is 

Zf < zt~ 5.07. (78) 

In the right panel, for (32 = 4.02, we obtain Zf < zt — 3.47. At the two-sigma lower limit, P2 = 1.55, the 
acceleration extends to Zf ~ 8. 

By contrast, it is easy to see that in dark-energy models of more conventional type, i.e. uncoupled to dark 
matter, with frozen equation of state and matter density f2m, the acceleration starts at the redshift 

zacc = [{Sw^ - 2)(n„ - l)/n„]^''(='-«"'*' - 1. (79) 

This value, for all < 2/3 (i.e. for a present accelerated regime), is always smaller than unity if ilm = 
0.3±0.1. Therefore, a (future) unambiguous measurement of the expansion rate at 2: > 1 could be a powerful 
method to distinguish between coupled and uncoupled models of dark energy. 

The plots of Fig. 1 refer to our "phenomenological" set of parameters, and in particular to the duration of 
the freezing epoch (possibly relevant to the solution of the coincidence problem). The as constraint provides 
however interesting information also on the set of "fundamental" parameters of the string effective action 
(0), used for our model of dilatonic dark energy. 

By eliminating Ci/c2 — -y/2/3/i2, and fixing qo, i.e. f)2, as before, we can plot indeed the as constraint in 
the plane spanned by the variables and Vq. The result is shown in Fig. 2, again for (32 = 2.35 (left panel) 
and for P2 = 4.02 (right panel). We have restored the required Planck length factors, and given the potential 
Vo in units of eV*. The allowed region is within the upper and lower white curves, as before. The lower 
bound on derived from (3i, at fixed (32, is satisfied for all the range of values illustrated in the picture. 

We note that the values of parameters used in |^] for a particular numerical integration of the string 
cosmology equations are well compatible with the above bounds. It is also important to stress that, since 
the value of is naturally expected in the range 10^-10^ §, the allowed mass scale of the dilaton potential, 
Vq^'^, turns out to be fixed in a rather narrow region around (10^^-10^^) eV, even in the case of an early 
start of the freezing epoch. This result confirms that, as already pointed out in |^, a realistic dark-energy 
scenario seems to require a rather high degree of accuracy in determining the scale of the dilaton potential. 
This (fine-tuning?) aspect of the potential is however a common problem of all scalar-field models of dark 
energy. 

It should be mentioned, to conclude this section, that we have performed an additional observational test 
of our class of dilatonic dark-energy models by comparing with the COBE data the slope of the low-multipole 
Ce spectrum induced by the SW and ISW effect. Using a simple Gaussian likelihood distribution we have 
concluded that, at the confidence level of two sigma, the predicted slopes are compatible with the data for 
all the allowed region of parameter space, so that no (significant) additional constraints are generated. 



V. THE FARTHEST SUPERNOVA 



As a last check on the viability of the present dilatonic dark-energy model we have considered the con- 
straints imposed by the most distant Type la supernova |Q known so far, SN1997ff, for which a very recent 
assessment of lensing magnification has increased the apparent magnitude by 0.34±0.12 mag This leads 
to a final distance modulus (i.e. to a difference of apparent and absolute magnitude) oim — M = 45.49±0.34 
mag. 
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Figure 1: Curves at constant ug (irom eq. (|76DJ, and /?2 fixed at the values 2.35 (top) and 4.02 (bottom). 
The allowed region is below the dashed line, and within the upper and lower white curves, corresponding to 
the limiting values of eq. ([t^). 



Using the definition adopted in |l2| ] , this result can be expressed in the following way. The luminosity 
distance is 

= + (80) 

from which one defines 

A(m - M) = 5[logio(rfL(^)) - logio(dLo(2))l, (81) 
where dLo(z) = z(z + 2)/(2_ffo) is the luminosity distance for Milne's model, i.e. an hyperbolic empty 




1 1.5 2 2.5 3 3.5 4 
log(c2) 



Figure 2: Curves at constant as, and (32 Sxed at the values 2.35 (top) and 4.02 (bottom). The allowed region 
is below the straight line representing the consistency bound Zf < zt, and within the upper and lower white 
curves corresponding to the limiting values of eq. (|7?[). 
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universe {flm = f^A = 0). For the supernova SN1997ff one obtains A(m - M) ^ -0.15 ± 0.34 at .z ~ 1.755, 
in good agreement with a ACDM model characterized today by Q,m = 0.35, SIa ~ 0.65 jl3[. For such a 
model, the Universe aX z — 1.755 is already well within the decelerated epoch, which starts around z = 0.548 
(see eq. |7|). 

This does not imply, however, that all models which are accelerated at large z are ruled out (even without 
mentioning the still unclear experimental uncertainties of such a supernova detection). Let us calculate 
indeed the luminosity-distance along a stationary Ptj,, for a spatially flat geometry. Neglecting 

for the moment the baryon contribution, the Friedmann equation is 

= Hi [n^ia/aa)-''^ + ^^(a/ao)-^"'] = Hl{a/ao)-^^, 

where, in particular, w = W2 (see eq. for our freezing regime. The corresponding luminosity-distance 
(for w ^ 2/3) is: 

dz' 2{l + z) 

For an accelerated evolution with w < 2/3 (i.e., in our case, go > 1), and for large z, we have (11 ~ 2^"^™''^, 
while for Milne's cosmology cIl ~ z^ . It follows that, at large z and for any w > 0, the Milne model 
always provides larger distances (and thus larger apparent magnitudes) than a model of stationary evolution. 
As a consequence, a negative value of A(m — M), referred to Milne, does not necessarily corresponds to 
deceleration. 

For a more precise illustration of this important point we have plotted in Fig. 3 the distance modulus 
A(m — M) for the accelerated freezing phase of our dilatonic dark-energy model. We have numerically 
integrated the luminosity-distance functions, including baryons, for the two particular values P2 = 2.35, f32 ~ 
4.02 already used in the previous flgures (it may be useful to recall that, when f^c is flxed to 0.3, these two 
values of (32 correspond to W2 — 0.5 and W2 — 0.4, respectively). It can be seen from the picture that, in 
both cases, the curves representing the cosmic evolution of our model, although deeply inside the accelerating 
regime, remain well within one sigma from the (lensing-corrected) SN1997ff data, while providing, at the 
same time, a reasonable fit of the binned data of all the other supernovae. 



(1 + ^) 



(82) 



VI. CONCLUSION 



In this paper we have considered a phenomenological model of dark-energy-dark-matter interactions based 
on the infinite bare-coupling limit of the superstring effective action. The dilaton, rolling down an exponen- 
tially suppressed potential, plays the role of the cosmic field responsible for the observed acceleration, and 
drives the Universe towards a final configuration dominated by a comparable amount of kinetic, potential 
and CDM energy density. 

The effective dilatonic coupling to dark matter switches on at late enough times (i.e., large enough bare 
coupling), and affects in a significant way the post-equivalence cosmological evolution. The time-dilution 
of the dark-matter density, in particular, is first slightly enhanced (during the dragging phase) and then 
considerably damped (during the freezing phase) with respect to the standard decay law. The large- 
angle fluctuation scales relevant to the observed CMB anisotropies reenter the horizon during the dragging 
epoch, and exit the horizon again during the freezing epoch. In spite of this unconventional evolution, the 
growth of the matter-density perturbations may be large enough to match consistently present observations. 

The predicted value of the (smoothed out) density contrast crs , compared with data obtained from cluster 
abundance, deflnes a signiflcant allowed region in the parameter space of the given class of dilatonic dark- 
energy models. The analysis of such an allowed region provides two main results. 

The flrst is that the bounds on the past-time extension of the accelerated (freezing) epoch are signiflcantly 
weaker than in conventional dark-energy models (uncoupled to dark matter, with frozen equation of state). 
The establishment of the freezing regime, in our class of dilatonic models, is allowed long before the present 
epoch (up to z ~ 5), thus providing (in principle) a further relaxation of the coincidence problem, by 
extending the present cosmological configuration not only in the far future, but also towards the past. 

The possibility of very early {z > 1) accelerated evolution is indeed a typical signature of such a class of 
dilatonic models, useful in principle to discriminate it from other (uncoupled) dark energy models, hopefully 
on the grounds of future observational data. It is important to stress, to this respect, that the farthest type 
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Figure 3: The distance modulus m—M, referred to Milne's cosmology, for a ACDM model (dashed curve) and 
for dilatonic dark-energy model during the freezing epoch, for two values of P2 (thick curves when baryons 
are included, thin curves without baryons). The points at z < 1 represent the binned data from all the 
high-redshift la supernovae known so far. The datum plotted at z — 1.755 represents the lensing-corrected 
SN1997S . All data are derived from ^ 



la supernova so far observed is at z ~ 1.7, and is perfectly compatible with an accelerated Universe already 
at that epoch, provided the data of the magnitude-redshift diagram are consistently fitted by the accelerated 
kinematics of dilatonic models. 

The second results concerns the parameters of the (non-perturbative) dilaton potential appearing in the 
strong-bare-coupling regime of the string effective action. The dilaton mass scale Vb, for an efficient and 
realistic dark energy scenario, appears in such a context to be tightly anchored to a value very near to the 
present Hubble curvature scale. A small deviation of Vo from the required value is enough to remove the 
predictions of the dilatonic model from the region of parameter space allowed by the erg data. 

This means that, under the assumption that the dilatonic models discussed in this paper provide the 
correct explanation of the observed cosmic acceleration, the measurements of the density contrast erg, be- 
sides their obvious astrophysical importance, would also acquire an interesting high-energy significance for 
providing an indirect (parameter-dependent) measurements of the dilaton mass scale. 

We note, finally, that astrophysical observations may provide several additional constraints on dilatonic 
dark-energy models. For instance, the clustering evolution of sources at high redshifts may constrain directly 



17 



the freezing growth exponent m2; as already mentioned, the baryon bias that develops during freezing is 
also observable, at least in principle [^; finally, further constraints can be derived from a computation of 
the full multipole spectrum of the CMB radiation (see e.g. I ^if for a recent study of the CMB constraints on 
coupled dark-energy models with power-law potentials). Preliminary results seem to confirm the conclusions 
of this paper, but a detailed discussion of these new constraints is postponed to a future work. 
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